The effect of harmonized emissions on aerosol properties in global models – an AeroCom experiment

The effects of unified aerosol sources on global aerosol fields simulated by different models are examined in this paper. We compare results from two AeroCom experiments, one with different (ExpA) and one with unified emissions, injection heights, and particle sizes at the source (ExpB). Surprisingly, harmonization of aerosol sources has only a small impact on the simulated diversity for aerosol burden, and consequently optical properties, as the results are largely controlled by model-specific transport, removal, chemistry (leading to the formation of secondary aerosols) and parameterizations of aerosol microphysics (e.g. the split between deposition pathways) and to a lesser extent on the spatial and temporal distributions of the (precursor) emissions. The burdens of black carbon and especially sea salt become more coherent in ExpB only, because the large ExpA diversity for these two species was caused by few outliers. The experiment also indicated that despite prescribing emission fluxes and size distributions, ambiguities in the implementation in individual models can lead to substantial differences. These results indicate the need for a better understanding of aerosol life cycles at process level (including spatial dispersal and interaction with meteorological parameters) in order to obtain more reliable results from global aerosol simulations. This is particularly important as such model results are used to assess the consequences of specific air pollution abatement strategies

An AeroCom initial assessment – optical properties in aerosol component modules of global models

The AeroCom exercise diagnoses multi-component aerosol modules in global modeling. In an initial assessment simulated global distributions for mass and mid-visible aerosol optical thickness (aot) were compared among 20 different modules. Model diversity was also explored in the context of previous comparisons. For the component combined aot general agreement has improved for the annual global mean. At 0.11 to 0.14, simulated aot values are at the lower end of global averages suggested by remote sensing from ground (AERONET ca. 0.135) and space (satellite composite ca. 0.15). More detailed comparisons, however, reveal that larger differences in regional distribution and significant differences in compositional mixture remain. Of particular concern are large model diversities for contributions by dust and carbonaceous aerosol, because they lead to significant uncertainty in aerosol absorption (aab). Since aot and aab, both, influence the aerosol impact on the radiative energy-balance, the aerosol (direct) forcing uncertainty in modeling is larger than differences in aot might suggest. New diagnostic approaches are proposed to trace model differences in terms of aerosol processing and transport: These include the prescription of common input (e.g. amount, size and injection of aerosol component emissions) and the use of observational capabilities from ground (e.g. measurements networks) or space (e.g. correlations between aerosol and clouds)

Analysis and quantification of the diversities of aerosol life cycles within AeroCom

Simulation results of global aerosol models have been assembled in the framework of the AeroCom intercomparison exercise. In this paper, we analyze the life cycles of dust, sea salt, sulfate, black carbon and particulate organic matter as simulated by sixteen global aerosol models. The diversities among the models for the sources and sinks, burdens, particle sizes, water uptakes, and spatial dispersals have been established. These diversities have large consequences for the calculated radiative forcing and the aerosol concentrations at the surface. The AeroCom all-models-average emissions are dominated by the mass of sea salt (SS), followed by dust (DU), sulfate (SO_4), particulate organic matter (POM), and finally black carbon (BC). Interactive parameterizations of the emissions and contrasting particles sizes of SS and DU lead generally to higher diversities of these species, and for total aerosol. The lower diversity of the emissions of the fine aerosols, BC, POM, and SO_4, is due to the use of similar emission inventories, and does therefore not necessarily indicate a better understanding of their sources. The diversity of SO_4-sources is mainly caused by the disagreement on depositional loss of precursor gases and on chemical production. The diversities of the emissions are passed on to the burdens, but the latter are also strongly affected by the model-specific treatments of transport and aerosol processes. The burdens of dry masses decrease from largest to smallest: DU, SS, SO_4, POM, and BC. The all-models-average residence time is shortest for SS with about half a day, followed by S_O4 and DU with four days, and POM and BC with six and seven days, respectively. The wet deposition rate is controlled by the solubility and increases from DU, BC, POM to SO_4 and SS. It is the dominant sink for SO_4, BC, and POM, and contributes about one third to the total removal rate coefficients of SS and DU species. For SS and DU we find high diversities for the removal rate coefficients and deposition pathways. Models do neither agree on the split between wet and dry deposition, nor on that between sedimentation and turbulent dry Deposition. We diagnose an extremely high diversity for the uptake of ambient water vapor that influences the particle size and thus the sink rate coefficients. Furthermore, we find little agreement among the model results for the partitioning of wet removal into scavenging by convective and stratiform rain. Large differences exist for aerosol dispersal both in the vertical and in the horizontal direction. In some models, a minimum of total aerosol concentration is simulated at the surface. Aerosol dispersal is most pronounced for SO4 and BC and lowest for SS. Diversities are higher for meridional than for vertical dispersal, they are similar for a given species and highest for SS and DU. For these two components we do not find a correlation between vertical and meridional aerosol dispersal. In addition the degree of dispersals of SS and DU is not related to their residence times. SO_4, BC, and POM, however, show increased meridional dispersal in models with larger vertical dispersal, and dispersal is larger for longer simulated residence times

Midlatitude atmospheric circulation responses under 1.5C and 2.0C warming and implications for regional impacts

This study investigates the global response of the midlatitude atmospheric circulation to 1.5◦C and 5 2.0◦C of warming using the HAPPI “Half a degree Additional warming, Prognosis and Projected Im- pacts” ensemble, with a focus on the winter season. Characterizing and understanding this response is critical for accurately assessing the near-term regional impacts of climate change and the benefits of limiting warming to 1.5◦C above pre-industrial levels, as advocated by the Paris Agreement of the United Nations Framework Convention on Climate Change (UNFCCC). The HAPPI experimental 10 design allows an assessment of uncertainty in the circulation response due to model dependence and internal variability. Internal variability is found to dominate the multi-model mean response of the jet streams, storm tracks and stationary waves across most of the midlatitudes; larger signals in these features are mostly consistent with those seen in more strongly forced warming scenarios. Signals that emerge in the 1.5◦C experiment are a weakening of storm activity over North America, an inland 15 shift of the North American stationary ridge, an equatorward shift of the North Pacific jet exit, and an equatorward intensification of the South Pacific jet. Signals that emerge under an additional 0.5◦C of warming include a poleward shift of the North Atlantic jet exit, an eastward extension of the North Atlantic storm track, and an intensification on the flanks of the Southern Hemisphere storm track. Case studies explore the implications of these circulation responses for precipitation impacts in the 20 Mediterranean, western Europe and on the North American west coast, paying particular attention to possible outcomes at the tails of the response distributions. For example, the projected weakening of the Mediterranean storm track emerges in the 2◦C warmer world, with exceptionally dry decades becoming five times more likely

Black carbon vertical profiles strongly affect its radiative forcing uncertainty

The impact of black carbon (BC) aerosols on the global radiation balance is not well constrained. Here twelve global aerosol models are used to show that at least 20% of the present uncertainty in modeled BC direct radiative forcing (RF) is due to diversity in the simulated vertical profile of BC mass. Results are from phases 1 and 2 of the global aerosol model intercomparison project (AeroCom). Additionally, a significant fraction of the variability is shown to come from high altitudes, as, globally, more than 40% of the total BC RF is exerted above 5 km. BC emission regions and areas with transported BC are found to have differing characteristics. These insights into the importance of the vertical profile of BC lead us to suggest that observational studies are needed to better characterize the global distribution of BC, including in the upper troposphere

Aerosols at the poles: an AeroCom Phase II multi-model evaluation

Atmospheric aerosols from anthropogenic and natural sources reach the polar regions through long-range transport and affect the local radiation balance. Such transport is, however, poorly constrained in present-day global climate models, and few multi-model evaluations of polar anthropogenic aerosol radiative forcing exist. Here we compare the aerosol optical depth (AOD) at 550 nm from simulations with 16 global aerosol models from the AeroCom Phase II model intercomparison project with available observations at both poles. We show that the annual mean multi-model median is representative of the observations in Arctic, but that the intermodel spread is large. We also document the geographical distribution and seasonal cycle of the AOD for the individual aerosol species: black carbon (BC) from fossil fuel and biomass burning, sulfate, organic aerosols (OAs), dust, and sea-salt. For a subset of models that represent nitrate and secondary organic aerosols (SOAs), we document the role of these aerosols at high latitudes. The seasonal dependence of natural and anthropogenic aerosols differs with natural aerosols peaking in winter (sea-salt) and spring (dust), whereas AOD from anthropogenic aerosols peaks in late spring and summer. The models produce a median annual mean AOD of 0.07 in the Arctic (defined here as north of 60° N). The models also predict a noteworthy aerosol transport to the Antarctic (south of 70° S) with a resulting AOD varying between 0.01 and 0.02. The models have estimated the shortwave anthropogenic radiative forcing contributions to the direct aerosol effect (DAE) associated with BC and OA from fossil fuel and biofuel (FF), sulfate, SOAs, nitrate, and biomass burning from BC and OA emissions combined. The Arctic modelled annual mean DAE is slightly negative (−0.12 W m−2), dominated by a positive BC FF DAE in spring and a negative sulfate DAE in summer. The Antarctic DAE is governed by BC FF. We perform sensitivity experiments with one of the AeroCom models (GISS modelE) to investigate how regional emissions of BC and sulfate and the lifetime of BC influence the Arctic and Antarctic AOD. A doubling of emissions in eastern Asia results in a 33 % increase in Arctic AOD of BC. A doubling of the BC lifetime results in a 39 % increase in Arctic AOD of BC. However, these radical changes still fall within the AeroCom model range

An AeroCom assessment of black carbon in Arctic snow and sea ice

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng g−1 for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng g−1 for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng g−1. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model–measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60–90° N) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07–0.25) W m−2 and 0.18 (0.06–0.28) W m−2 in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W m−2 for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic

Evaluation of Black Carbon Estimations in Global Aerosol Models

We evaluate black carbon (BC) model predictions from the AeroCom model intercomparison project by considering the diversity among year 2000 model simulations and comparing model predictions with available measurements. These 5 model-measurement intercomparisons include BC surface and aircraft concentrations, aerosol absorption optical depth (AAOD) from AERONET and Ozone Monitoring Instrument (OMI) retrievals and BC column estimations based on AERONET. In regions other than Asia, most models are biased high compared to surface concentration measurements. However compared with (column) AAOD or BC burden retreivals, the models 10 are generally biased low. The average ratio of model to retrieved AAOD is less than 0.7 in South American and 0.6 in African biomass burning regions; both of these regions lack surface concentration measurements. In Asia the average model to observed ratio is 0.6 for AAOD and 0.5 for BC surface concentrations. Compared with aircraft measurements over the Americas at latitudes between 0 and 50 N, the average model is a 15 factor of 10 larger than observed, and most models exceed the measured BC standard deviation in the mid to upper troposphere. At higher latitudes the average model to aircraft BC is 0.6 and underestimates the observed BC loading in the lower and middle troposphere associated with springtime Arctic haze. Low model bias for AAOD but overestimation of surface and upper atmospheric BC concentrations at lower latitudes 20 suggests that most models are underestimating BC absorption and should improve estimates for refractive index, particle size, and optical effects of BC coating. Retrieval uncertainties and/or differences with model diagnostic treatment may also contribute to the model-measurement disparity. Largest AeroCom model diversity occurred in northern Eurasia and the remote Arctic, regions influenced by anthropogenic sources. 25 Changing emissions, aging, removal, or optical properties within a single model generated a smaller change in model predictions than the range represented by the full set of AeroCom models. Upper tropospheric concentrations of BC mass from the aircraft measurements are suggested to provide a unique new benchmark to test scavenging and vertical dispersion of BC in global models.JRC.H.2-Climate chang

Aerosol indirect effects

Aerosol indirect effects continue to constitute one of the most important uncertainties for anthropogenic climate perturbations. Within the international AEROCOM initiative, the representation of aerosol-cloud-radiation interactions in ten different general circulation models (GCMs) is evaluated using three satellite datasets. The focus is on stratiform liquid water clouds since most GCMs do not include ice nucleation effects, and none of the model explicitly parameterises aerosol effects on convective clouds. We compute statistical relationships between aerosol optical depth (tau a) and various cloud and radiation quantities in a manner that is consistent between the models and the satellite data. cloud droplet number concentration (N d) compares relatively well to the satellite data at least over the ocean. The relationship between (tau a) and liquid water path is simulated much too strongly by the models. This suggests that the implementation of the second aerosol indirect effect mainly in terms of an autoconversion parameterisation has to be revisited in the GCMs. A positive relationship between total cloud fraction (fcld) and tau a as found in the satellite data is simulated by the majority of the models, albeit less strongly than that in the satellite data in most of them. In a discussion of the hypotheses proposed in the literature to explain the satellite-derived strong fcld–tau a relationship, our results indicate that none can be identified as a unique explanation. Relationships similar to the ones found in satellite data between tau a and cloud top temperature or outgoing long-wave radiation (OLR) are simulated by only a few GCMs. The GCMs that simulate a negative OLR - tau a relationship show a strong positive correlation between tau a and fcld. The short-wave total aerosol radiative forcing as simulated by the GCMs is strongly influenced by the simulated anthropogenic fraction of tau a, and parameterisation assumptions such as a lower bound on Nd. Nevertheless, the strengths of the statistical relationships are good predictors for the aerosol forcings in the models. An estimate of the total short-wave aerosol forcing inferred from the combination of these predictors for the modelled forcings with the satellite-derived statistical relationships yields a global annual mean value of −1.5±0.5Wm−2. In an alternative approach, the radiative flux perturbation due to anthropogenic aerosols can be broken down into a component over the cloud-free portion of the globe (approximately the aerosol direct effect) and a component over the cloudy portion of the globe (approximately the aerosol indirect effect). An estimate obtained by scaling these simulated clearand cloudy-sky forcings with estimates of anthropogenic tau a and satellite-retrieved Nd–tau a regression slopes, respectively, yields a global, annual-mean aerosol direct effect estimate of −0.4±0.2Wm−2 and a cloudy-sky (aerosol indirect effect) estimate of −0.7±0.5Wm−2, with a total estimate of −1.2±0.4Wm−2

Evaluating the climate and air quality impacts of short-lived pollutants

This paper presents a summary of the work done within the European Union's Seventh Framework Programme project ECLIPSE (Evaluating the Climate and Air Quality Impacts of Short-Lived Pollutants). ECLIPSE had a unique systematic concept for designing a realistic and effective mitigation scenario for short-lived climate pollutants (SLCPs; methane, aerosols and ozone, and their precursor species) and quantifying its climate and air quality impacts, and this paper presents the results in the context of this overarching strategy. The first step in ECLIPSE was to create a new emission inventory based on current legislation (CLE) for the recent past and until 2050. Substantial progress compared to previous work was made by including previously unaccounted types of sources such as flaring of gas associated with oil production, and wick lamps. These emission data were used for present-day reference simulations with four advanced Earth system models (ESMs) and six chemistry transport models (CTMs). The model simulations were compared with a variety of ground-based and satellite observational data sets from Asia, Europe and the Arctic. It was found that the models still underestimate the measured seasonality of aerosols in the Arctic but to a lesser extent than in previous studies. Problems likely related to the emissions were identified for northern Russia and India, in particular. To estimate the climate impacts of SLCPs, ECLIPSE followed two paths of research: the first path calculated radiative forcing (RF) values for a large matrix of SLCP species emissions, for different seasons and regions independently. Based on these RF calculations, the Global Temperature change Potential metric for a time horizon of 20 years (GTP20) was calculated for each SLCP emission type. This climate metric was then used in an integrated assessment model to identify all emission mitigation measures with a beneficial air quality and short-term (20-year) climate impact. These measures together defined a SLCP mitigation (MIT) scenario. Compared to CLE, the MIT scenario would reduce global methane (CH4) and black carbon (BC) emissions by about 50 and 80 %, respectively. For CH4, measures on shale gas production, waste management and coal mines were most important. For non-CH4 SLCPs, elimination of high-emitting vehicles and wick lamps, as well as reducing emissions from gas flaring, coal and biomass stoves, agricultural waste, solvents and diesel engines were most important. These measures lead to large reductions in calculated surface concentrations of ozone and particulate matter. We estimate that in the EU, the loss of statistical life expectancy due to air pollution was 7.5 months in 2010, which will be reduced to 5.2 months by 2030 in the CLE scenario. The MIT scenario would reduce this value by another 0.9 to 4.3 months. Substantially larger reductions due to the mitigation are found for China (1.8 months) and India (11–12 months). The climate metrics cannot fully quantify the climate response. Therefore, a second research path was taken. Transient climate ensemble simulations with the four ESMs were run for the CLE and MIT scenarios, to determine the climate impacts of the mitigation. In these simulations, the CLE scenario resulted in a surface temperature increase of 0.70 ± 0.14 K between the years 2006 and 2050. For the decade 2041–2050, the warming was reduced by 0.22 ± 0.07 K in the MIT scenario, and this result was in almost exact agreement with the response calculated based on the emission metrics (reduced warming of 0.22 ± 0.09 K). The metrics calculations suggest that non-CH4 SLCPs contribute ~ 22 % to this response and CH4 78 %. This could not be fully confirmed by the transient simulations, which attributed about 90 % of the temperature response to CH4 reductions. Attribution of the observed temperature response to non-CH4 SLCP emission reductions and BC specifically is hampered in the transient simulations by small forcing and co-emitted species of the emission basket chosen. Nevertheless, an important conclusion is that our mitigation basket as a whole would lead to clear benefits for both air quality and climate. The climate response from BC reductions in our study is smaller than reported previously, possibly because our study is one of the first to use fully coupled climate models, where unforced variability and sea ice responses cause relatively strong temperature fluctuations that may counteract (and, thus, mask) the impacts of small emission reductions. The temperature responses to the mitigation were generally stronger over the continents than over the oceans, and with a warming reduction of 0.44 K (0.39–0.49) K the largest over the Arctic. Our calculations suggest particularly beneficial climate responses in southern Europe, where surface warming was reduced by about 0.3 K and precipitation rates were increased by about 15 (6–21) mm yr−1 (more than 4 % of total precipitation) from spring to autumn. Thus, the mitigation could help to alleviate expected future drought and water shortages in the Mediterranean area. We also report other important results of the ECLIPSE project